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ABSTRACT 

It has recently been shown that galaxy formation models within the ACDM cosmology 
predict that, compared to the observed population, small galaxies (with stellar masses < 10^^ 
Mq) form too early, are too passive since z ~ 3 and host too old stellar populations at z = 0. 
We then expect an overproduction of small galaxies at z > 4 that should be visible as an excess 
of faint Lyman-break galaxies. To check whether this excess is present, we use the MORGANA 
galaxy formation model and GRASIL spectro-photometric + radiative transfer code to gener- 
ate mock catalogues of deep fields observed with HST-ACS. We add observational noise and 
the effect of Lyman-alpha emission, and perform color-color selections to identify Lyman- 
break galaxies. The resulting mock candidates have plausible properties that closely resemble 
those of observed galaxies. We are able to reproduce the evolution of the bright tail of the 
luminosity function of Lyman-break galaxies (with a possible underestimate of the number of 
the brightest i-dropouts), but uncertainties and degeneracies in dust absorption parameters do 
not allow to give strong constraints to the model. Besides, our model shows a clear excess with 
respect to observations of faint Lyman-break galaxies, especially of zgso ^ 27 F-dropouts at 
z ^ 5. We quantify the properties of these "excess" galaxies and discuss the implications: 
these galaxies are hosted in dark matter halos with circular velocities in excess of 100 km 
s~^, and their suppression may require a deep re-thinking of stellar feedback processes taking 
place in galaxy formation. 
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1 INTRODUCTION 

Despite the large observational and theoretical efforts of the last 
decade, the formation of galaxies in the context of the ACDM 
cosmology is still an open problem. In fact, the wealth of multi- 
wavelength data on deep fields, coupled with wide local surveys, 
has produced new constraints on number densities, stellar masses, 
star formation rates and metallicities of galaxies in the redshift 
range from to ~ 7. From the phenomenological point of view, 
a remarkable behaviour of observed galaxies is their "downsizing" 
trend in star formation rates, stellar masses, chemical enrichment, 
reconstructed age of stellar populations, reconstructed duration of 
the m ain star formation episode (see references in Fontanot et^ 
For "downsizing" here we mean that more massive galaxies 
tend to be older than less massive ones in terms of either stellar 
ages or mass assembly, a behaviour that is at variance with that of 
dark matter halos. 



These constraints are a strong ch allenge to theoretical m odels 
of galaxy formation. In a recent paper, [pontanot et^ ( |2009|) com- 
pared the predictions of thr e e galaxy formation mod els l lWang et"ail 
l2008l : [Monaco etai]l2007l : ISomerville e7al]|2008h with observa- 
tional data which showed downsizing trends. This analysis allowed 
the authors to highlight a significant discrepancy between mod- 
els and data that pertains to relatively small galaxies, with stellar 
masses in the range 10^ — 10^"'^ M©. In all three models these 
galaxies are formed at high redshift and are rather old and passive 
at late times, while the population of real galaxies in the same stel- 
lar mass range builds up at low z, growing in number by a factor 
of ~ 3 from 2: ~ 1 to the present time, and shows higher specific 
star formation rates up to 2: ~ 3 — 4. Similar conclusions are ob- 
tained considering the age of stellar populations in local galaxies: 
real galaxies host much younger stars on average than model galax- 
ies. The consistency of results coming from stellar masses, star for- 
mation rates and reconstructed stellar ages makes this trend robust. 
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Previous comparisons, focused on stellar mass downsizing, showed 
that this discrepancy is not limited to the three models tested in that 
paper but extends to many semi-analytic and N-bo dy codes (see 
the comparisons to models in iFontana et al. 2006: Cirasuolo et alj 
bOOSLlMarchesim et alfeOOSl) . 

This is a severe discrepancy because, as we will discuss later 
in this paper, it cannot be easily solved by a simple fine-tuning of 
parameters. Also, it is not directly related to the well-known ten- 
dency of semi-analytic models to produce steep luminosity func- 
tions, because all the model stellar mass functions have roughly 
the correct slope at 2: = 0, or to the well-known problem of 
overquenching of satellite galaxies, which ar e too red and pas - 
sive in models (see, e.g., the recent analysis o f iKimm et al.ll2008b . 
Clearly the solution of this problem amounts to understanding why 
do small galaxies form so early and why do they stop forming stars 
later The second problem is most likey a byproduct of the first one: 
once too many small galaxies are produced at high redshift, strong 
stellar feedback must be invoked to suppress their star formation at 
lower redshift, so as to recover the correct number density at 2; = 0. 
Then, the crucial objects that are at the origin of this discrepancy 
are the small star-forming galaxies at 2; > 3. Unfortunately, with 
presently available facilities only the rest-frame UV radiation of 
such galaxies is accessible; these are visible as distant faint Lyman- 
break galaxies. Their name takes origin from the technique used to 
select them, also known as "dropout technique". This technique se- 
lects galaxies over a specific redshift interval by detecting the drop 
of the UV flux observed in correspondence with the Lyman limit at 
A < 912 A, the wavelength below which the ground state of neutral 
hydrogen may be ionized. 

In this paper we use the MORGANA -1- GRASIL model of galaxy 
formation to generate predictions of B-, V- and i-dropouts to be 
compared with the re sults of (i) the GOODS -S p hotometric and 
spectr oscopic data by iGiavalisco et al.l ( l2004t) and IVanzella et al.l 
( |2008|) . (ii) the GOODS-MUSIC sample tarazian et al."2006') with 
photometric redshifts, (iii) the results of Bouwens et al. ( 2003), 
where many HST-ACS deep fields were analyse d and scaled to 
the d epth of the Hubble Ultra Deep Field (HUDF. Iseckwith et al.l 
I2OO6I) . 

Several previous theoretical studies addressed the proper- 
ties of Lyman-break galaxies using both sem i-analytical models 
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2008: Lee et al. 2008: Guo & White 2008) and cosmological simu- 
latio ns (Nagamine 2002: Weinberg et al. 2002; Harford & Gn edinI 
2003i:.Nagamine et al 2004; .Finlator et al..2006 ; Night et al..26o6F 



Nagamine et al. I l2008h . These papers checked how models within 



the ACDM cosmology are able to reproduce the properties of 
Lyman-break galaxies and their clustering, often reporting satis- 
factory agreement . Desp i te some author s (namely IPoli et alfcOOll 
I2OO3I : iNight et al.l I2OO6I ; iFinlator et all l2006h noticed that mod- 
els tend to overproduce faint high-redshift Lyman Break Galax- 
ies with respect to observations, no quantitative claim had been 
possible at that time, given the uncertainties on the observational 
constraints. As an example. IFinlator et al. I ( l2006h found that their 
predicted steep faint-ends of the luminosity function were still in 
agreement with the available data, which spanned a range of val- 
ues from —2.2 to —1.6. T his wide range has recently been revised 
by teouwens et al.ll2007l) . so that the most recent determinations 
give a slope ~ 1.7 ± 0.1 from z ~ 2 to z ~ 6, which implies a 
disagreement with the theoretical predictions. 

The aim of this paper is a careful revision of this compari- 
son, in the light of the improved observational constraints, based 



on the generation of reliable and improved mock galaxy catalogues 
through state-of-the-art semi-analytical and spectro-photometric 
models. This approach will allow us to test the impact of the selec- 
tion criteria, as well as different theoretical ingredients (with par- 
ticular emphasis on the modeling of dust attenuation). Moreover, 
we aim to include this comparison in a cosmological context, and 
to discuss its implications for cuiTent models of ga laxy formation 
and the so-called "downsizing" scenario (see e.g. iFontanot et al.l 
(ilooi)). 

The plan of the paper is the following. Section|2]presents the 
model, which contains some improvements with respect to previ- 
ous versions. Section [3] describes the adopted procedure to select 
Lyman-break galaxies from the simulated galaxy catalogues and to 
tune the free parameters connected to dust extinction. Section |4] 
presents the comparison of model and data in terms of number 
counts, redshift and color distributions and luminosity functions 
of B,V and i-dropouts, while Section [5] quantifies the properties 
of the model "excess" galaxies and discusses future observational 
tests that can help shedding light on this population. Section|6]gives 
the conclusions. 



2 THE MORGANA -1- GRASIL MODEL 
2.1 MORGANA 

The MOdel for the Rise of GAlaxies aNd Agn s (MORGANA) is de- 
scribed in full detail in [Monaco et al.l ( 120071) . T he model foll ows 
the general scheme of semi-analytic models (see lBaugl]|2006l . for 
a review): here we only mention how the main processes are im- 
plemented, (i) Each dark matter halo hosts a central galaxy and 
many satellite galaxies, contained within dark matter substruc- 
ture that has not yet been destroyed by orbital decay due to dy- 
namical friction. Baryons are subdivided into stellar halo, bulge 
and disc components. Each component has cold, hot and stellar 
phases, (ii) Merge r trees of dark matter halos are obtained using the 
PINO CCHIO tool llMonaco et al.l 2002; [Monaco, Theuns & Taffonil 
2002; iTaffoni et al. 2002). (iii) After a merging of dark matter 



halos, galaxy merg ing times are computed using the results of 
iTaffoni et alj ( l2003h . which are tuned on N-body simulations that 
include dynamical friction, tidal stripping and tidal shocks, (iv) 
At each merger a fraction of the satellite stellar mass is moved to 
the diffuse stellar co mponent of the stellar halo, as explained in 
[Monaco et al. I ( l2006h . (iv) The evolution of the baryonic compo- 
nents and phases is performed by numerically integrating a system 
of equations for all the mass, energy and metal flows, (v) The inter- 
galactic medium infalling on a dark matter halo is shock-heated to 
a specific thermal energy equal to 1.2 times that of dark matter, (vi) 
At each time-step of the numerical integration, the hot gas com- 
ponent is assumed to be in hydrostatic equilibrium with the dark 
matter halo, (vii) Cooling flows are computed with a new cooling 
model, tested against simulations of isolated halos in Iviola et al.l 
( i2008h . which takes into account mass and energy input from feed- 
back sources, (viii) The cooling gas is let infall on the central 
galaxy on a dynamical time-scale. It is divided between disc and 
bulge according to the fraction of the disc that lies within the half- 
mass radius of the bulge, (ix) The gas infalling on the disc keeps 
its angular momentum; disc sizes are computed with an extension 
of the Mo, Mao & White (1998) model that includes the contribu- 
tion of the bulge to the disc rotation curve, (x) Disc instabilities and 
major mergers of galaxies lead to the formation of bulges. In minor 
mergers the satellite mass is given to the bulge component of the 
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larger galaxy, (xi) Star formation and f eedback in discs are inserted 
fol lowing the model o f lMonacol J2OO4I) . while in bulge components 
the lKennicuti ill 9981) law is used. In both cases, hot gas is ejected to 
the halo (in a hot galactic wind) at a rate equal to the star-formation 
rate, (xii) Feedback from star formation heats the hot halo phase. 
When its temperature gets significantly higher than the virial one, 
it leaves the dark matter halo in a galactic hot super-wind. The 
matter is re-acquired later, when the virial temperature of the de- 
scendant halo gets as high as that of the hot superwind gas. (xiii) 
In star-forming bulges cold gas is ejected in a cold galactic wind 
by kinetic feedback due to the predicted high level of turbulence 
driven by SNe. Analogously to the hot superwind, cold halo gas can 
be ejected out of the dark matter halo if its kinetic energy is high 
enough, (xiv) Accretion of gas onto massive black holes (starting 
from small seeds present in all galaxies) is connected to the abil- 
ity of cold gas to loose angular momentum by some (unspecified) 
mech anisms driven by star f ormation. This is explained in full de- 
tail in lFontanot et al.l 1 I2OO6I) . (xv) The kind of feedback from the 
AGN depends on the accretion rate in Eddington units: whenever 
this is hig her than I per cent the AGN can trigger a massive galaxy 
wind, (see [Monaco & Fontanolll2005h which leads to the complete 
removal of Inter-Stellar Medium (ISM) from the galaxy, while at 
lower accretion rates the energy is ejected through jets that feed 
back on the hot halo gas; this can lead to quenching of the cooling 
flow. 

2.2 The improvements with respect to the earlier version 

The version of the mo del used here present s a few improvements 
with respect to that of [Monaco et alj ( |2007|) . The most important 
difference lies in the treatment of star formation in the bulge com- 
ponent. In the earlier version, the size of the starburst, used to com- 
put e the gas surface density and then the star formation rate through 
the lKennicuttI ( Il998l) relation, was assumed to be equal to the bulge 
effective radius Rb- As a consequence, a starburst declines very 
slowly because the residual gas has a decreasingly low surface den- 
sity, and then a long star formation time-scale. This is a clearly 
naive estimate, because the dissipational gas settles in the centre of 
the galaxy. In this version we assume that the size of the starburst 
is determined by the level of turbulence that the gas acquires. The 
velocity dispersion (JcoM due to turbulence is: 
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where the parameter ctq is set to 60 km s ^ . iFontanot et al. 

show that this value allows a good reproduction of the evolution of 

the AGN population. The time-scale of star formation is: 
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where Ecom is the (cold) gas surface density. The size of the 
startburst is determined as the radius for which the turbulence- 
driven velocity dispersion equals that of the host bulge computed 
at the same radius. For the Young mass profile used in the model, 
the bulge velocity curve is roughly constant (assuming a value 
~ 0.707Vbuigc, where Vbuigc is the circular velocity of the bulge) 
down to 0.04 times the effective radius Rb, where it quickly drops 
to zero. Consequently, the size of starbursts in bulges is estimated 
as: 
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where (JcoidiRb) is computed using Rb to determine the surface 
density (equation |2](. The radius -Rgas is not allowed to be smaller 
than OMRb. 

A seco nd difference is that we optimized this version of the 
model for a lChabrierl j2003h mass function instead of the Salpeter 
one used before. To achieve the same level of agreement with 
data as with Sal peter we perfo r m the following parameter changes 
(please refer to iMonaco et alj ( |2007|) for all details): (i) the stel- 
lar mass per supernova is M*,sn = 84 Mq, the restoration rate 
is /rest ~ 0.44, while the metal yield per generation is set to 
Y = 0.016; (ii) the preferred values of some parameters change 

to: Tlquench = 0.3, Tldyn = 0.3, /wind = 2.0, /th,B = 0.1, 

/th,D = 0. 32, /iet.o = 3.0. Qu asar-triggered galaxy winds are 
present as in lFontanot et alj ( [2006h ; dark matter halo concentrations 
are not rescaled to improve the fit of the baryonic TuUy-Fisher re- 
lation. Moreover, the threshold for disc instability is relaxed from 
eiimit ~ 0.9 to 0.7; this allows to remove the excess of small 
bulges present in the previous version of the model (see Figure 7 of 
iMonaco et"al]|2007l) . 

Finally, to slow down the merger-driven evolution of the most 
massive galaxies, at each merger a fraction /scatter = 0.8 of the 
stellar mass of each satellite is given to the halo diffuse stellar 
component; howeve r, this is done only at z < 1, as suggested by 
IConrovetal] ( l2007l) . This is not a physicall y acceptable solution of 
the problem raised bv lMonaco et al.l ( |2006|) . but has the advantage 
of giving realistic results. 

2.3 GRASIL 

For each galaxy modeled by MORGANA the corrisponding UV-to- 
radio SED is compute d with the spectr o-photometric + radiative 
transfer code GRASIL jSilva et alj|l99^ please refer to this paper 
for a detailed description of this code). According to this model, 
stars and dust are distribuited in a bulge (King profile) + disc (radial 
and vertical exponential profiles) axisymmetric geometry. Three 
different dusty environments are taken into account: (1) dust in 
interstellar HI clouds heated by the general interstellar radiation 
field of the galaxy (the "cirrus" component), (2) dust associated 
with star-forming molecular clouds and HII regions, and (3) cir- 
cumstellar dust shells produced by the windy final stages of stellar 
evolution. These environments have different importance in differ- 
ent galactic systems at various evolutionary stages. Therefore the 
residual gas fraction in the galaxy is divided into two phases: the 
dense giant molecular clouds (MCs) and the diffuse medium ("cir- 
rus"). The fraction of gas in the molecular phase, fmoi, is a model 
parameter. The molecular gas is subdivided into clouds of given 
mass and radius: the stars are assumed to be born within the op- 
tically thick MCs and gradually to escape from them, as they get 
older, on a time-scale tcsc- To compute the emitted spectrum of 
the star-forming molecular clouds a radiative transfer code is used. 
The diffuse dust emission is derived by describing the galaxy as 
an axially symmetric system, in which the local dust emissivity is 
consistently computed as a function of the local field intensity due 
to the stellar component. 

For each galaxy modelled by MORGANA, star formation his- 
tories, together with gas and metal contents at the "observation" 
time for both bulge and disc component, are then given as an in- 
put to GRASIL. Unobs cured stellar SEDs are generated by using 
agam a IChabrieil ( I2OO3I) IMF with stellar masses in the limit [0.15- 
120] Mq. In the radiative transfer computation the effect of dust 
is limited to absorption; dust emission, which is very expensive 
from the point of view of computational time, is irrelevant in the 
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rest-frame UV spectral regions sampled here. Absorption from the 
Inter-Galactic Medium (IGM) i s taken into account by applying to 
the redshifted SEDs t he lMadaul l ll995h attenuation, as explained in 
iFontanot et al]j2007h . 

Here below the GRASIL parameters that are not provided 
by MORGANA and which play an important role in the dust 
p arametrization are brie fly discussed. We also recall the values used 
in lFontanot et al.l ( |2007|) to compare with if -band and sub-mm se- 
lected galaxies. 

(i) tesc represents the tim e-scale for young O-B stars to es- 
cape from their parent MCs. In lFontanot et alj ( IiOOTI) we fixed the 
value of tesc to lO'yr.. Thi s is an intermediate value between those 
found bv lSilva et al.l l ll998l) to well describe the SED of spirals (~ 
a few Myr) and starburst (~ a few 10 Myr), and it is of the order 
of the estimated destruction time-scale of MCs by massive stars, 
(ii) The cold gas mass provided by MO RGANA is s ubdivided be- 
tween the dense and diffuse phases. In lFontanot et a l. (2007) we 
fixed the fraction of gas in MCs /moi to O.5.. (iii) The total (i.e., 
cirrus -I- MCs) dust content of a galaxy ISM is obtained by the resid- 
ual gas mass and the dust-to-gas mass ratio 5 which is assumed to 
scale linearly with the metallicity {5 = 0.45Z). The optical depth 
of MCs depends on their ma ss and radius o n ly thr ough the rela- 
tion r oc SMmc /'r\ic In iFontanot et al.l ( |2007|) we assumed 
Mmc = 10^ Mq, tmc = 16pc. (iv) The bulge and disc scale 
radii for stars and gas are given by MORGANA, while the disc 
sc ale-heights, h*i and hi for stars and dust, respectively, are set as 
in lFontanot et al.l ( l2007h to 0.1 the corresponding scale radii. The 
dust grain size distribution and abundances are set so as to match 
the me an Galactic extinction curve and emissivity as in lSilva et al.l 
( Il998l) and lVega et "Zl | |2005|) and are not varied here. 

Because the rest-frame UV emission is especially sensitive to 
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high redshift). In this paper we will explore the effect of changing 
the values of the two parameters tesc and fmoi keeping fixed the 
other ones. These parameters, in fact, by definition, are those which 
can play a fundamental role in the attenuation of the UV emission 
from young massive stars. This will be discussed in more detail in 
the next section. 



2.4 Lyman-a emission line 

Many Lyman-break galaxies show a prominent Lyman-a emis- 
sion line, which may influence galaxy colors. This is especially 
true at the highest redshift, where the i-z color can get signifi- 
cantly bluer, so that the selection of a galaxy as an i-dropout is 
shifted to higher redshift. Modeling the effect of such emission 
is very difficult due to the resonant scattering of Lyman-a pho- 
tons by neutral hydrogen which implies that even a small amount 
of dust is sufficient to absorb many emitted photons and attenuate 
the Lyman-a line. The observed lines are generally ascribed to the 
presence of galactic winds which allow Lyman-a photons to escape 
the galaxy after a few scatterings. Radiative transfer calculations of 
Lyman-a through winds have shown that thi s process can explain 
the observed asymmetric line profiles (e. g . IShaDlevetal]|2003l ; 
IVanzella et al ] |2006l ; IVerhamme et"al]|2008l) . Due to this complex- 
ity, we prefer to make an empirical estimate of the possible influ- 
ence of the Lyman-a line on galaxy colors and number counts. 

To this aim, we add a Lyman-a emission line to the produced 
spectra, with specified EW . Given the relatively coarse wavelength 
sampling used (sufficient to compute magnitudes over standard fil- 



ters), this line is added to one single wavelength bin. With the 
VLT/F0RS2 spectroscopic campaign of the GOODS-S field the 
number of Lyman-break candidates with observed Lyman-a emis- 
sion lines is 25 over 85 for the B-dropouts, 19 over 52 for the V- 
dropouts, 24 over 65 for the i-dropouts. The media n EW for the 
three cases are 15li| 221^^ A and 641^°" A jVanzella et al] 
I2OO9I) . These numbers should be taken with a grain of salt, because 
most dropout candidates (mainly the fainter ones) are safely con- 
firmed only when a Lyman-a emission line is found, so it is not 
easy either to compute an unbiased estimate of the fraction of emit- 
ters with respect to confirmed dropouts or to assess the complete- 
ness limit on the EW distribution. As a rough guess we assume that 
1 /3 of the galaxies have emission lines, and assume for the EW 
either the median values reported above for the three dropout cat- 
egories, applying them at redshifts [3, 4.5], [4.5, 5.5] and [5.5, 00], 
or the y-dropout value of 22 A for all e mitters. The se values are 
roughly consistent with those assumed bv lBouwens e t al. (2007) in 
their analysis. 

We will show result of this procedure only for one combina- 
tion of GRASIL parameters; analogous results are obtained in other 



3 SIMULATIONS OF DEEP FIELDS 

To run MORGANA we use merger trees from a 1 000^ particles 
PINOCCHIO run in a 200 Mpc box with the WMAP3 ( ISpergell26o^ 
cosmology Q.o = 0.24, f^A = 0.76, Q.b = 0.0456, Ho = 73 
km/s/Mpc, (jg = 0. 8, Us = 0.96 . The power spectrum was com- 
puted using the fit bv lEisenstein & Hu ( 1998). The particle mass in 
this run is 2.84 ■ 10* Mq, so the smallest used dark matter halo is 
1.42 • 10^" M0 (50 particles) and the smallest resolved progenitor 
is 2.84 • 10^ Mq (10 particles). 

To simulate a deep field, we need to pass from a time sampling 
of galaxies in the box to a pencil beam where time is translated to 
redshift, then compute galaxy SEDs in the observer frame. 

We transform the output from MORGANA into a pencil beam 
as follows. The model gives as an output, for each galaxy and in a 
fixed time grid, its physical properties like mass and metal content 
of each phase in each component, bulge and disc si zes and circular 
velocities, bulge and disc star formation rates (see IFontanot et al.l 
I2OO7I for details). The time grid length has been set to At = 10 
Myr, and the run has been stopped at z = 3. 

As specified in Paper I, the mass function of dark matter halos 
is sparse-sampled, computing the evolution of at most 300 halos 
per bin of 0.5 dex in log mass. This results in an oversampling of 
small satellite galaxies compared to the central galaxies of simi- 
lar stellar mass. To correct this oversampling, satellites are further 
sparse-sampled so as to be as abundant, at the final redshift (2 = 3), 
as the corresponding central galaxies of the same stellar mass. An 
approximately constant number of galaxies in logarithmic intervals 
of mass is thus obtained from this procedure. This way, the num- 
ber of produced galaxies is anyway very high, because each galaxy 
can be "observed" at each output time, i.e. each 10 Myr A further 
sparse sampling (1 over 100) is then applied to all galaxies, so as to 
select a sufficiently large sample of roughly 35,000 galaxies. Each 
of these sparse samplings define a weight for the galaxies, equal 
to the inverse of the selection probability. All the statistical quan- 
tities are then computed by weighte d sums over the galax ies. Our 
procedure, explai ned in full detail in Fontanot et alj ( |2007|) . is dif- 
ferent from, e.g., iKitzbichler & Whitd ( I2OO7I) . where information 
on the spatial position of the halos in the simulated box is retained. 
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Our procedure has the advantadge of being simpler and, by aver- 
aging over the box (redshift does not depend on the position of the 
galaxy in the box), giving results insensitive to large scale structure, 
while the procedure of iKitz bichler & White ( 2007) is necessary to 
address, for instance, the clustering of high-redshift galaxies. 

The galaxy sample is then transformed into a pencil beam 
by assigning to each galaxy "observed" at time t a redshift z and 
proper, luminosity and diameter distances r(z), dh and do- The 
comoving volume associated to each galaxy is computed as the co- 
moving area of the simulation box times the Ar(2:) distance asso- 
ciated to the cosmic time interval At of 10 Myr. The solid angle 
to be associated to the galaxy is that subtended by the box at the 
"observed" redshift z. 

The apparent (AB) magnitudes of the galaxies are computed 
by convolving the resulting redshifted SEDs from GRASIL with 
the ACS filters -B435, Veoe, ^775 and zsso (in the following for sim- 
plicity we drop the pedices from the magnitudes, with the exception 
of 2850 in order to avoid confusion with the redshift). 

In order to better reproduce the observed number counts of 
Lyman-break galaxies and their color distributions, photometric 
scatter has been associated to each galaxy magnitude. Model- 
ing this scatter is important because of the tendency for fainter, 
lower signal-to-noise sources to scatter into the selection through a 
Malmquist-like effect. To quantify this scatter it is necessary, first 
of all, to obtain an estimation of the flux error associated to the 
apparent magnitude of the observed galaxies. 

Observational errors are assigned to magnitudes on the 
basis of the ACS GOODS catalogue, version 2.0. Details 
of the ACS observations, as well a s major features of th e 
GOODS project, can be found in iGiavalisco et al. I ( l2004h : 
additional information about the latest v2.0 release of the 
GOODS ACS images and source catalogues can be found 
at |http : / /www .stsci.edu /science/goods /v2.0| and 
will be described in detail in an upcoming paper. Here we only 
mention that we use the best guess for the total magnitude for com- 
parison with the model, while galaxy colours of real galaxies are 
always computed using isophotal magnitudes. Source detection for 
this catalog has been performed on the z^zo band, the isophotes de- 
fined by the zsso image have been then used as apertures for all 
other bands. Using this catalogue we compute for each band {B, 
V, i, Z85o) the average signal-to-noise ratio as a function of the AB 
apparent magnitude of the observed galaxies. For B, V and i bands 
this has been done down to l-a detection limit, because it is im- 
portant to set proper upper limits in order to select faint drop-outs. 
For the reference band, 2;850, where by definition the sources are 
detected above 5-a, we use the i band error to have a good guess of 
the l-a error. Then, given the AB magnitude of a model galaxy, its 
flux / is perturbed with a Gaussian-distributed random error with 
width Gf equal to the corresponding average l-a error flux com- 
puted above from the GOODS catalogue. Fluxes are finally trans- 
formed back to magnitudes. Whenever the theoretical magnitude is 
larger (fainter) than that corresponding to S/N = 1, which means 
that such objects are not detectable in that band {S/N < 1), we 
put their magnitude equal to the corresponding l-a limit. The new 
magnitude represents therefore the lower limit, (upper limit if we 
are looking at the flux of the object), at l-a to the real value. The 
corresponding l-a limit magnitudes are 30.45 (-B435), 30.56 (Vaoa) 
and 30.19 (2775 and zsso). 

This procedure is tuned to reproduce ACS deep fields, so 
we will compare our predictions to the spectroscopic sa mple of 
GOODS-S Lyman-break galaxies jVanzeUa et"aLl l2009h, to the 
GOODS-MUSIC sample with photometric redshifts 1 Grazian et al.l 



|2006|) . which is based on the version 1.0 of the ACS GOODS cata- 
logue including deep J and K data of VUT-ISAAC and mid-IR data 
of Spitzer-I RAC, with photometric redshifts, and with the results of 
iBouwensTt al. ( 2007). The latter authors measured number counts 
and luminosity functions of Lyman-break galaxies in various ACS 
deep fields, homoge nizing the samples and scaling all results to the 
depth of the HUDF l lBeckwith et"al]|2006l) . 

Selection of Lyman break galaxies is then applied to the appar- 
ent magnitudes computed as explained above. We use the following 
criteria respectively for the B-, V- and i-dropout ( IGiavalisco et al.l 
l2004h : 

{B~V > 1.1) AND [B-V > {V- Z85o) + 1.1] 

AND {V - zsso) < 1.6 (4) 

{[V - i > 1.5 + 0.9{i - zs5o)] OR{V-i> 2)} 

AND {V-i> 1.2) AND (i - ^850) < 1.3 (5) 

(i - z85o) > 1.3 AND S/N{B) < 1 AND S/N{V) < 1 (6) 



3.1 TXining grasil 

The prediction of galaxy apparent magnitudes critically depends on 
dust attenuation, for which several uncertain parameters must be 
specified in GRASIL. Among all the parameters, we identify two of 
them as the crucial ones, namely /moi, the fraction of gas in molec- 
ular clouds, and tcsc, the time after which OB stars get ou t of the 
highly obscured molecular cloud. In iFontanot et alj ( l2007h values 
of /moi ~ 0.5 (typical of the Milky Way) and tosc = 10 Myr were 
used to make predictions in the K and sub-mm bands. However, 
when applied to our (much more dust-sensitive) Lyman-break can- 
didates they give dust attenuations of more than 2 mag, definitely 
larger, e.g., than the v alues of 1.4, 1.15 and 0.8 mag assumed by 
iBouwens et al.l (l2007i) for B-, V- and i-dropouts. Such parameter 
values may not be suitable for compact, gas-rich starbursting galax- 
ies as our Lyman-break candidates are. In fact, higher molecular 
fractions are o bserved in galaxies where the ISM is more pressur- 
ized (see, e.g.. lBlitz & Rosolowskvll2004h : this would imply much 
lower attenuation by diffuse dust in the ISM (cirrus component). 

We then run grasil using all combinations of /moi ~ 0.5, 0.90, 
0.95, 1.0 and tosc = 10, 3, 1 Myr. Tabled gives a fist of all the 
runs performed. We use the model described in Section ETI with 
quasar-triggered galaxy winds (std), but perform also a few runs 
without winds (nowind) to check the effect of this (very uncertain) 
ingredient of the model. For sake of compl eteness we also show 
res ults for the same model (old) presented in lMonaco et al. I( l2007t) 
and lFontanot et alj ( [2007h . with GRASIL run with a Salpeter IMF, 
but with /moi = 0.50 and tcsc ~ 3 Myr, a combination which gives 
a best-fit of the Lyman-break luminosity functions. 

For one of the models (std.f095.e03) we perform two 
runs with the addition of the Lyman-a emission line, either as- 
suming a redshift-dependent EW (std . f 95 . e03 . lyal) or 
assuming a constant EW (std.f095.e03.lya2), as motivated 
in Sectionl2!4l 



4 RESULTS 

Figure [T] shows the color-color plots for the selection of the model 
B-, V- and i-dropouts, according to the criteria defined above. Here 
and in the next figure we use the std.f095.e03 model; here 
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Figure 1. Color-color diagrams of B-, V- and i-dropouts for the std . f 95 . e03 model. (Blue) filled squares are galaxies with redshift in the specific 
interval of the selection, that is, [3.5-4.5], [4.5-5.5] and [5.5-6.5] for B-, V- and i-dropouts respectively. The open (cyan) squares represent the interlopers (see 
text for details), while the (red) open triangles ai'e upper-limits in B,V and i. For sake of clarity only the central galaxies (see Section lTTl are here considered, 
(satellite galaxies have very similar colors). 
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Figure 2. Selection function for B-. V- and i-dropouts for the std.f095.e03 model. The fraction of selected galaxies is computed in bins of redshift and 
apparent magnitude, and is reported as a tone of gray, as indicated in the right of each panel. 



for sake of clarity only central galaxies with zgso < 28, the res- 
olution limit of the model, are used (satellite galaxies have very 
similar colors). All galaxies in the redshift intervals of the selec- 
tion [3.5-4.5], [4.5-5.5] and [5.5-6.5] for B-y- and i-dropouts re- 
spectively, are shown as full blue squares, or red triangles if they 
are upper limits. Thes e are the same redshift intervals identified by 
iBouwens et alj l l2007h as typical of Lyman-break galaxies at mean 
redshift of z ~ 4, z ~ 5 and 2 ~ 6. The few interlopers falling in 
the selection areas are marked as cyan open squares. For interlop- 
ers we mean galaxies with redshift in the interval [3.0-3.4] for the 
B-dropouts, [4.0-4.4] for the V^-dropouts and [5.0-5.4] for the i- 
dropouts. Clearly, most galaxies in the relevant redshift interval are 
selected. Those that scatter out of the selection region are mostly 
faint galaxies with upper limits in the band from which they drop 
out; their colors are made less red by the upper limits. Figure |2] 
shows the selection function as a function of redshift and i (B- 
dropouts) or 2350 (V^- and i-dropouts) magnitude. Selected galaxies 
clearly tend to reside in the expected redshift range, though some 
B- and V-dropouts are lost at magnitudes fainter than ~ 26. 

We first compare model predictions of B-, V- and i-dropouts 



with observed number counts, which are the directly observable 
quantity when complete spectroscopic samples are not available. 

Figure [3]shows the effect of changing the GRASIL parameters 
tesc and /moi compared to the standard model without dust attenua- 
tion (black dot-dashed lines). We consider both the effect of chang- 
ing fcsc keeping fixed the value of /moi (upper panel), and the ef- 
fect of varying /moi at fixed tcsc (bottom panel). Model predictions 
for unabsorbed spectra drop at a magnitude ~ 27.5 (fainter for the 
i-dropouts). This is the typical magnitude of the central galaxy con- 
tained in the smallest resolved dark matter halo (1.42 ■ 10^" Mq), 
fainter galaxies would be hosted in smaller, unresolved halos and 
are thus under-represented. It is worth noting that HUDF data are 
deeper than this limiQ Clearly this limit magnitude increases when 



^ It would be easy too reach deeper magnitudes by running the model on 
a smaller box. However, here we are interested in understanding whether 
the same model that produces small galaxies that are too old at z = 
)Fontanot et"ai]l2009l) is also over-producing faint star-forming galaxies at 
high redshift. To achieve this, the requirement is that observations probe at 
least as deep as the model's limiting magnitude. 
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Figure 3. Number counts of B-, V- and t-dropouts. (Red) points are from Bouwens et al. (2007). Upper panels: models with fi-^ol = 0.95, std.f095.el0 
(cyan solid lines), std.f095.e03 (green long-dashed Hues) and std.f095.e01 (blue dashed lines) compared to the model without dust attenuation (black 
dot-dashed lines). Lower panels: models with fesc = 0.3, std.f050.e03 (cyan solid lines), std.f090.e03 (blue dashed lines) and std.f095.e03 
(green long-dashed lines) and std.fl00.e03 (magenta dotted lines) also compared to the model without dust attenuation (black dot-dashed lines). 



dust attenuation is considered. When compared to Bouwens et al.'s 
observations (open red points in the Figure), number counts in ab- 
sence of dust attenuation clearly overshoot the data (by about 2 
mag). This shows how a correct modeling of attenuation is funda- 
mental to best-fit the data. This feature clearly limits the predictive 
power of the model, because small variations of dust parameters 
within the allowed range give very large differences in the results. 

The figure shows that the two parameters are degenerate; we 
identify the two combinations (tesc, /moi) = (3 Myr, 0.95) (model 
std.f095.e03) and (1 Myr, 0.90) (model std. f 090 . eOl) 
as the best-fitting ones. This degeneration concerns the extinction, 
we in fact, expect to find differences in the IR due to the different 
distributions of temperature and optical depth of the two phases. 
The attenuations given by these best-fit models are roughly simi- 
lar to the ones used by Bouwens et al. (1.4, 1.15 and 0.8 mag for 
B-, V- and i-dropouts) but are systematica lly higher for brighte r 
galaxies. This is in line with the findings of IShapley et al.| 120011): 
these authors reported for a large sample of Lyman-break galaxies 
at 2 ~ 3 an average E{B — V) attenuation of 0.15 that correlates 
with magnitude, brighter objects being more absorbed. Expressed 
in terms of the same E{B — V) quantity, our attenuation is 0.1 at 
the bright end and declines to ~ 0.05 for the faint objects. 

Figure |4] shows number counts for the two best-fit models. As 
expected, while the bright end is well reproduced with its redshift 



evolution (with a modest excess for the B-dropouts and possibly a 
dearth of bright i-dropouts), a clear excess of faint objects is present 
at Z850 27, especially for the ^/-dropouts. 

The inclusion of the Lyman-a line does not change much the 
predicted number counts of Lyman-break galaxies, with the excep- 
tions of the i-dropouts, which are especially sensitive to the line. 
Figure |5] shows the effect of including the Lyman-a line, as ex- 
plained in Section l2!4l to the std.f095.e03 best-fit model. As 
expected, the effect of the Lyman-a emission line is to make the 
1-^850 color bluer, thus shifting the entrance in the selection crite- 
rion to slightly higher redshift and decreasing the number counts. 
This decrease is significant at the bright end, which is then under- 
estimated by the model (in one case it also removes the kink at 
bright luminosity, which is due to the presence of few, high-weight 
galaxies - see Section[3}, while the faint end is not strongly influ- 
enced. This conclusion does not change much with the choice of 
the EW. The inclusion of the Lyman-a emission line performed 
here is clearly tentative, but this test shows that our results are ro- 
bust with respect to the inclusion of Lyman-a emission line under 
reasonable assumptions. 

We then concentrate on the std.f095.e03 model, to un- 
derstand to what extent model Lyman-break galaxies resemble the 
real ones. 

Figure|6]shows the B-V, V-i and i-zgso colours of model and 
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Table 1. Runs performed for this paper, with values of GRASIL parameters 
(see Section l23] for a brief description of the parameters) and indication of 
the MORGANA model used. 



real galaxies as a function of redshift. Observed galaxies are shown 
both from the GOODS-S catalogue, with spectroscopic redshift, 
and for the GOODS-MUSIC catalogue, with photometric redshift, 
which clearly has many more galaxies. Only the galaxies brighter 
than 2850 < 26.2 are here considered, because of the 90 per cent 
completeness limit of the GOODS-MUSIC catalogue. As already 
noticed in Fig. [T] model galaxies present colours that closely re- 
semble the observed ones, though they tend to be slightly redder, 
especially in B-V at z~4 (at least compared to GOODS-MUSIC). 
Analogous conclusions can be drawn by considering the other best- 
fit model, std.f090.e01. This means that colours cannot be 
used to constrain the parameters. This conclusion is strenghtened 
by noting that the strong reddening of all colors with redshift is 
driven by IGM absorption. 

Figure |7] shows the redshift distributions of model Lyman- 
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Figure 5. Number counts of «-dropouts. (Red) points are from Bouwens et 
al. (2007), the (green) solid line gives the predictions of the best-fit model 
std.f095.e03, the (black) dashed and the (cyan) long-dashed lines cor- 
respond to the same combination of parameters but with the contribution of 
Lya emission line in the two different prescription discussed in Section lT4l 
(black: std.f095.e03. lyal.cyan: std. f 95 . e03 . Iya2). 



break galaxies. We compare these predictions with the spectro- 
scopic redshift distribution of GOODS-S gal axies (red solid lines) 
which constitute an incomplete sample (see IVanzella et al" 
[2009;, for details) and with the photometric redshift distributions of 
GOODS-MUSIC (black dashed lines) and consider only galaxies 
brighter than ^sso = 26.2, the range where the GOODS-MUSIC 
sample is complete. The agreement is good though not perfect; this 
can be seen by comparing the average and standard deviation of 
the predicted redshift distributions, (3.71 ± 0.35; 4.88 ± 0.27; 
5.67 ± 0.22 ; for B-,V- and i^dropout respectively), with those 
derived from lGiavalisco et al.l ( [2004b . (3.78 ± 0.34; 4.92 ± 0.33; 
5.74 ± 0.36; for the three candidates respectively), the GOODS- 
S spectroscopic sample (3.78 ± 0.37; 4.98 ± 0.39; 5.92 ± 0.20) 
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Figure 6. B-V, V-i and i-^gso colors of galaxies with zg5o < 26.2 as a function of redshift. Left panels: GOODS-S galaxies with spectroscopic redshift. 
Mid panels: GOODS-MUSIC galaxies with photometric redshifts. Black arrows indicate that the magnitude estimation is based on an upper limit at l-cr. Right 
panels: model galaxies (std . f 95 . e03). 



and GOODS-MUSIC photometric sample (3.81 ± 0.40; 4.95 ± 
0.37; 5.74 ± 0.10); however sample variance in the 150 arcmin^ 
GOODS-S field is not negligible, so the level of agreement is con- 
sidered satisfactory, especially if we consider that a confirmation of 
the validity of the redshift distribution is represented by the consis- 
tency between number counts and luminosity functions. Anyway, 
since the uncertainty related to the redshift measurements increases 
at higher redshifts where the detection is more difficult, this com- 
parison should be taken with grain of salt. In particular, the number 
of available photometric redshifts of i-dropout galaxies is very low, 



because only few such galaxies are bright enough to be detected in 
J and K bands. 

ISouwens et alj l l2007l) estimated the luminosity function of 
Lyman-break galaxies by scaling all observations to the HUDF (5- 
a) depth of 29 in the zgso band. Assuming that the B-dropout tech- 
nique selects all star-forming (and not severely obscured) galaxies 
in the redshift interval 3.6 < z < 4.5, they obtained selection 
functions of V- and i-dropouts by degrading images of low-z B- 
dropout galaxies to reproduce observational biases at higher red- 
shift, then estimating how many of them would be selected. As 



10 Lo Faro et al. 



1.5 



Z 1 



0.5 



/ \ 
f-V 



1.5 



0.5 



J I I I I I I I I 

3 3.5 4 4.5 5 
z 









A 


\ / 


\ 
\ 






V 


' \ 
'--\. 

\ 
;i 




w 



4 4.5 5 
z 



5.5 6 



2 - 



I 1 



5.5 6 6.5 7 
z 
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Figure 8. Luminosity function of B-, V- and i-dropouts. (Red) points are the estimate by Bouwens et al. (2007), (green) solid and (blue) dashed lines give 
predictions for the std.f095.e03 and std.f090.e01 models. 



shown in figures[T]and|2l most model galaxies in that redshift range 
with high enough jzgso flux are B-dropouts, so it is fair to compare 
Bouwens et al.'s estimate to the luminosity function of all model 
galaxies in the redshift range 3.5 < z < 4.5, 4.5 < z < 5.5 
and 5.5 < z < 6.5. Figure[8]shows t he luminosity function of the 
two best-fit models, compared to the iBouwens et alj ( |2007|) deter- 
mination, adjusted to take into account the slightly different cos- 
mology. The agreement is again good at high luminosities (with a 
possible underestimate of the i-dropout bright tail), while a clear 
excess is visible in the faint tail, especially in the V^-dropouts. We 
stress that this result is different from the one of figure |4] relative 
to number counts, because the two comparisons follow two differ- 
ent approaches: in the first case the model tries to reproduce ex- 
actly the observed quantity, in the second case the observers try to 
reconstruct a quantity more directly comparable with models, i.e. 
the luminosity function. The good agreement of data and model in 
terms of both number counts and luminosity functions cofirms that 
model galaxies closely resemble the observed population and their 
selection function. Moreover, the model follows very well the pos- 
itive evolution (with time) of the knee of the luminosity functions, 
and this allows us to support Bouwens et al.' conclusion that we 



are indeed witnessing the hierarchical build-up of the population of 
young galaxies. 

Summing up, the MORGANA model reproduces the main ob- 
servational properties (number counts, colours, selection functions, 
redshift distributions and estimated luminosity functions) of bright 
Lyman-break star-forming galaxies at z > 3.5, (though a possible 
underestimate of the number of bright z-dropouts may be present); 
fainter galaxies, with zgso > 27, have acceptable properties but are 
over-produced. However, dust attenuation heavily depends on un- 
certain parameters of the ISM like tcsc and /moi, so that multiple 
solutions are equally acceptable. 



5 A POPULATION OF EXCESS GALAXIES 

Degeneracies in model predictions do not allow us to draw strong 
conclusions from the good agreement shown above between model 
and real Lyman-break galaxies at the bright end of the luminosity 
function. Indeed, parameter degeneracy is not due only to GRASIL; 
two MORGANA parameters, namely Jidyn and riqucnch (see sec- 
tion |2.2| l influence the star-formation rate density at z > 3 but give 
small difference at lower redshift, where more constraints are avail- 
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able. To strengthen this point, we show in Figure |9] the prediction 
for the luminosity function of a model run without quasar-triggered 
galaxy winds (nowind . f 9 5 . e03) and the old. f 090 . e03 
used in iFontanot et al. I l l2007h , where a Salpeter IMF is used in 
place of Chabrier, and a slightly different cosmology with as = 0.9 
is adopted. The last model is run on a 512'^ box with slightly worse 
mass resolution. The fit is very similar to our best fits in all cases, 
and this confirms that bright Lyman-break galaxies cannot give 
strong constraints to models as long as dust attenuation is so un- 
certain. The reason why quasar winds do not affect this observable 
is the following. In our imp lemented scheme for quasar-triggered 
winds jFontanot et alj200a) radiation from an accreting black hole 
can evaporate some cold gas in the star-forming bulge; if the puta- 
tive evaporation rate is higher than the star formation rate, then all 
the cold mass in the host bulge is removed. Quasar winds are then 
able to quench star formation in bulges (mergers), and this could 
lead to a suppression of star formation. However, due to the (mod- 
eled) delay between star formation and black-hole accretion, quasar 
winds take place when the star formation episode that has triggered 
accretion is almost over. The independence of the Lyman-break lu- 
minosity function on winds is a further demonstration of this. 

Previous papers frontana et al.' '2006|; ICirasuolo et alj 1200 Sl ; 
iMarchesini et alj200 8 ; Fontanot et al. 200^ have shown that MOR- 
GANA presents a deficit of massive galaxies at z > 3, when com- 
pared to the stellar mass function or to the K-band luminosity func- 
tion. At the same time, MORGANA is able to reproduce the sub- 
mm counts of z ~ 2 starburst galaxies jFontanot et al.|[2007h . and 
even shows a modest excess of bright B-dropouts (fig.|4]l. The two 
evidences of low stellar masses and a correctly high SFR are not 
contradictory: compared with similar models, MORGANA produces 
stronger cooling flows and then stronger bursts of star formation 
that are readily quenched by stellar feedback, so that the high lev- 
els of star formation are maintained for short times and do not pro- 
duce higher amounts of stellar mass. Moreover, the apparent un- 
derestimate of stellar masses by MORGANA m ay be at least in par t 
due to systematics in stellar mass estimates. iTonini et alj | |2008|) 
showed that inclusion of TP-AGBs in the Simple Stellar Popula- 
tion libraries give a significant downward revision of stellar masses 
at high redshifts or, equivalently, a boost to the j-iT-band m agnitude 
of model galaxies. According to lMarchesini et al] ( |2008|) . our stel- 
lar mass function at z ~ 3 remains low at the ~ 2a level, but 
when we compare stellar masses of our model Lyman-break galax- 



ies with estimates of IStark et al. I ( l2009l) . which take into account 
TP-AGBs, we find a slight over-prediction, in contrast with the evi- 
dences mentioned above. It is then difficult at this stage to assess to 
what extent a good match of the bright end of the luminosity func- 
tion of Lyman-break galaxies at 2 ~ 4 — 5 is in contrast with the 
apparent underestimate of the high end of the stellar mass function 
at 2 ~ 3. 

Once the bright tail is matched, the model shows a clear excess 
of faint dropouts. Is this excess robust? From the o ne hand, as ex- 
plaine d in the Introduction, the excess is expected: iFontanot et al.l 
demonstrated that galaxies with stellar masses (at 2 = 0) 
in the lO'"* — 10^*^ M© range are formed too early and are too 
passive since 2 ~ 3, and this is the main cause of the lack of 
(stellar mass and archaeological) downsizing in models. It follows 
that an excess of small star-forming galaxies must be present at 
higher redshift, and because the data go deeper than the complete- 
ness limit of the model, this excess must be visible. On the other 
hand, number counts at such deep magnitudes may well be incom- 
plete. The use of the HUDF, and the careful analysis of Bouwens et 
al. to determine the selection functions, makes it unlikely that the 
model excess is purely due to sample incompleteness. Moreover, 
the reconstructed faint slope of the luminosity function is ~ —1.7, 
in agreement with many other determinations both at 2 ^ 3 
jYoshida et al."2006':'Beckwith et al.|[2006l : lOesch et alj|2007l) and 
at 2 ~ 2 feeddy & Steidel JOO^, while the slope of model lumi- 
nosity functions is significantly steeper, in excess of —2.0. 

Is this excess an artifact of dust extinction? in other words, is it 
possible that MORGANA is underestimating the star formation rate 
of bright galaxies (as well as stellar masses) but our ad-hoc tun- 
ing of dust parameters absorbs this discrepancy, creating an excess 
of faint galaxies? Figure [Tol shows the E{B — V) attenuation as 
a function of UV rest-frame luminosity Misuo.ab (left panel) and 
gas-metallicity (right panel) for model V-dropout galaxies at 2 ~ 5 
and with Log{Mi,) > 7.75; central galaxies are highlighted as 
empty circles. The solid black lines represent the 'weighted' mean 
values of E{B — V) computed, respectively, for each bin of mag- 
nitude and gas-metallicity. For B- and i-dropouts we obtain a sim- 
ilar trend. The figure shows that the attenuation we obtain is only 
sli ghtly lower iE(B — V) ~ 0.1 in place of 0.15) than what found 
bv lShaplev et al.l 1 200 ih at lower redshift, (but we are not using the 
same SED template of those authors because every galaxy in our 
model has its own SED), while, as mentioned above, the attenua- 
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Figure 10. E{B — V) attenuation as a function of tlie UV rest-frame luminosity Migoo, AS (left panel) and gas-metallicity (referred to solar metallicity)(right 
panel). All model ^-dropout galaxies at redshift about 5 and with Log{Mi,) > 7.75, are here considered. Open circles refer to central galaxies. The black 
solid lines represent the 'weighted' mean values of E{B — V) computed, respectively, for each bin of magnitude and gas-metallicity. B- and i-dropouts show 
a similar trend. 



tio n in the observed i an d zgso bands is very similar to that used 
by iBouwens et al.l ( [2007h ; the excess remains if we attenuate our 
galaxies as suggested in that paper Moreover, the faint-end slope 
of model luminosity functions, in excess of —2.0, is inconsistent 
with the latest and very stable —1.7 value; with such a difference, 
an excess of faint galaxies will be obtained even if the bright tail 
is underestimated. This steep slope is in part due to to the fact that 
faint galaxies are less absorbed t han bright ones, soniething that is 
well known in literature (see, e.g., [Sanders & Mirabel 199a), and is 
consistent with the trend found by[Shaplev et al. (2001 ) at 2 ~ 3. 
We conclude that it is unlikely that this excess is due to a peculiar 
behaviour of dust attenuation as predicted by GRASIL. 

Other theoretical groups have compared their predictions of 
Lyman-break galaxies to data (see the reference list in the Intro- 
duction), but no group has clearly claimed to find a similar excess, 
with the exception of IPoli et al.l j200ll | 2003|), who comp ared the 
predictions of a semi-analytical model jMenci et al.l2002l) with the 
reconstructed luminosity function of I and /('-selected samples {not 
color selected). They reported an excess of faint / ~ 27 galaxies at 
lower reds hift, z ^ 3, which is ev en stronger than the one we find 
here. Also, iNight et alj ( |2006|) and lFinlator et al.l | |2006|) , extracting 
catalogues of Lyman-break galaxies from N-body simu lations, find 
rather steep luminosity functions, with slopes ~ —2. iNight et al.l 
( l2006l) discuss in some detail the implications of this steep slope, 
but claim broad consistency with available data; compared to the 
more recent determination of Bouwens et al., such steep slopes are 
nominally ruled out. There are two reasons why many groups may 
have not noticed such an excess: first, very deep fields like GOODS 
and the HUDF are required to reliably sample the zgso ~ 27 — 28 
magnitude range, and these have been made available only recently. 
Second, a simple (and over-simplistic) magnitude-independent at- 
tenuation makes the excess less visible, so a more advanced tool as 
GRASIL is needed to notice it. 

In the following of this section we try to assess whether the 
excess can be fixed simply by a better fine-tuning of model parame- 
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Figure 12. Stellar masses of the z = descendant galaxies of V-dropouts 
with 27 < 2850 < 28. 



ters, without spoiling the good match with local galaxies. To answer 
this question, we focus on the properties of the excess population. 
FigurefTTIshows stellar masses, star formation rates and dark matter 
circular velocities of B-, V- and i-dropout galaxies as a function of 
apparent zgso magnitude, for the best-fit model std.f095.e03; 
central galaxies are highlighted as empty circles. Excess galaxies 
with 2850 ~ 27 — 28 are characterized by star formation rates of 
order of < 10 M0 yr^'^, stellar masses growing from 10** - 10^ 
M0 at z ~ 6 to a factor of ~ 3 higher at z ~ 4, and the central 
ones are hosted in dark matter halos with masses in the range from 
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Figure 11. Properties of model Lyman-break galaxies as a function of 2350 magnitude, for the model std.f095.e03. Upper panels give stellar masses, 
mid panel give star formation rates, lower panels give dark matter halo circular velocities. Bigger open points refer to central galaxies. 



10^^ to lO^'^ M0, whose circular velocities are in excess of 100 
km s~^. Figure[T2] shows the stellar mass of the 2 = descendant 
galaxies of V-dropouts with 27 < 285 < 28. Roughly half of them 
do not grow much in mass since z ~ 5, thus falling in the stellar 
mass range of galaxies that are too old in models. The other half 
merge with larger galaxies to end u p at the knee of t he ste llar mass 
function. In many papers, like e.g. ISantini & et alj ( |200^ ). the ob- 
served specific star formation rates of massive galaxies is found to 
be higher than model predictions by a factor at least of ~ 4, which 
is apparently at variance with the good match of the growth of stel - 
lar mass in the same galaxies since z ~ 3 jpontana et al]|2006h . 



The puzzle may be solved if the excess population of small galax- 
ies contributes to inflate the growth of massive galaxies by minor 
mergers, thus compensating for the lower star formation rates. We 
propose this as a plausible explanation, although more stable esti- 
mates of star formation rates in galaxies are necessary to assess this 
problem. 

An important hint on the mechanism responsible for the ex- 
cess of faint galaxies comes from their metallicities Z. Obser- 
vations show that metallicities measured from emission lines de- 
cline with redshift, and this decline is st ronger for smal l er gala xies. 
This has been measured at z ~ 2 bv ISavaglio et and 
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Figure 13. Stellar mass - gas metallicity relation for model galaxies at 
z = 5. In this plot all galaxies in the box are used. Bigger open points 
refer to central galaxies. The two horizontal lines mai'k solar and 1/3 solar 
metallicities. We removed galaxies with less than 10 per cent gas fraction. 



by lErb et alj ( |2006|) . lErbl ( |2007|) . and at z ~ 3 by iMaiolino eTal 
(l2008n . but some measures exist at higher redshift jAndo et al 
20071) that s eem t o confirm the trend up to z ~ 6. Besides, 
Pettini et"al] ( 1200 ih report a yalue of 1/3 solar as the typical metal- 
licity of z ~ 3 Lyman-break galaxies. We show in figure [T3l the 
predicted metallicities of the gas component of galaxies at z = 5 
as a function of their galaxy stellar mass. We show gas metallicities 
because these are measured from emission lines. A mass-metal re- 
lation is yisible, but with a very high normalization, for which 10^*^ 
Mq galaxies are already enriched above 1 /3 solar, the typical gas 
metallicity of z ~ 3 Lyman-break galaxies. Moreover, a blob of 
small, highly enriched galaxies is visible, that is probably absent in 
the real world. These high metallicities justify the slightly redder 
colours of model galaxies and the shift toward lower redshift of the 
selection function, and may be at the origin of the underestimate of 
the number of bright i-dropouts. High metallicities at high redshift 
ar e a common problem of galaxy formation models, as pointed out 
in lMaiolino et al. | (l2008h: relat i vely h igh metallicities at z ~ 2 — 3 
are reported by iGuo & White (l2008l). and s tudyin g model predic- 
tions for the Milky Way lDe Lucia & Helmil fcoOSh noticed a dearth 
of low metallicity halo stars, due to excessive enrichment of small 
galaxies (before they become Milky Way satellite) at high red- 
shift. Can this mismatc h in gas metallicity cre ate the excess of faint 
galaxies? Observations jMaiolino et alj2008h suggest a downsizing 
trend for metal enrichment, with the stellar mass-gas metallicity re- 
lation steepening at high redshift. Our relation is rather flat, so that 
we are probably overestimating gas metallicity more in fainter ob- 
jects. If more metals give higher attenuation, this mismatch would 
lead to an artificiai flattening of the LF, thus decreasing the excess 
of faint Lyman-break galaxies rather than creating it. 

There is a broad consensus among modelers on the fact that 
AGN feedback is most effective in massive galaxies at low redshift, 
so the excess galaxies should be suppressed by stellar feedback. 
This should work in the direction of creating strong gas and metal 
outflows, so as to keep low both stellar mass and metallicity. Com- 



ing back to the question on whether this excess can be suppressed 
by fine-tuning model parameters, it is easy to do it by increasing the 
efficiency of stellar feedback, but at the cost of destroying the good 
fit of the stellar mass function at low redshift. Indeed, as clear in 
Figure[TT] such galaxies, though small, are hosted in relatively large 
dark matter halos with circular velocities well in excess of 100 km 
s~^. Because the efficiency of gas removal mainly sc ales with halo 
circul ar velocity (as found since the seminal work of lOekel & Silld 
Il986l) . this would be done at the cost of suppressing also star for- 
mation in those 10^^ Mq halos at z = in which star formation 
efficiency is highest, thus destroying the good agreement with the 
local stellar mass function (and luminosity functions as well). 

Stellar feedback is manily driven by SNe, and each SN pro- 
duces roug hly 10^^ erg of energy, which translates into 10 erg 
per Mq of stars formed if roughly 1 SN explodes each 100 Mq of 
solar masses. If a fraction e of SN energy is injected in the ISM 
instead of being quickly radiated away (e ~ 0.3 in our model), the 
velocity Weject at which a mass Mcjcct can be accelerated by a mass 
Mi, of stars formed is t'cjoct ~ 500 -^Z (e/0.3)A/*/Afcjcct km s^^. 
As a result, one Mq of stars can not eject much more than one Mq 
of gas out of a halo with circular velocity of ~ 200 km s^^ (and 
escape velocity higher by the canonical factor of \/2)- In order to 
limit the number and metal content of faint Lyman-break galaxies 
hosted in 100 km/s halos at z ~ 5, Uojoct must be higher than it 
is at z ~ 0, and this can be obtained as follows: (i) increasing en- 
ergy feedback efficiency e (if not the SN energy), (ii) increasing 
the number of SNe per solar mass of stars formed, (iii) decreasing 
halo circular velocities (or their number at fixed Vc). The first solu- 
tion would imply that feedback works differently in such starburst, 
with the paradoxical conclusion that feedback can be most efficient 
in the densest environments. Anyway, the gain factor is unlikely 
to be very high, because if 32 per cent of SN energy (in discs) is 
available for feedback, as assumed in this paper, the gain cannot 
be higher than a factor of three in energy (or Wcjcct)- Things would 
change if such SNe were more energetic than their low-redshift and 
high-metallicity counterparts. Solution (ii) would be obtained by 
assuming a different, more top-heavy IMF in these galaxies, while 
solution (iii) would require that dark matter is not completely cold 
and collisionless as assumed. 

In all cases, a solution of this discrepancy, and a proper re- 
production of stellar mass and archaeological downsizing, requires 
some deep change in the models, be it a more sophisticated feed- 
back scheme, more energetic SNe at high redshift, a varying IMF 
or a different dark matter theory. Moreover, this question is of great 
interest because, due to their steep luminosity function, such galax- 
ies would contain most produced metals, so they are very important 
contributors to metal pollution of the IGM. 

While theoretical progress can be achieved by exploring the 
ideas given above, observational breakthroughs will have to wait 
the next generation of telescopes. Indeed, while the photometric 
selection of Lyman-break galaxies is feasible with the current facil- 
ities up to magnitude zgso ~28, even though on a relatively small 
field like GOODS, the spectroscopic observations of relevant fea- 
tures like the break of the continuum blueward the Lya line and 
the detection of the ultraviolet absorption lines is still challenging 
for sources fainter than zgso ~ 25.5. Beyond this limit the main 
spectral feature is the Lyman a emission line (if present), often 
observed without continuum; only composite spectra allow to in- 
crease the signal-to-noise. Future observations in the next decade 
with ground-based 30-40m class telescopes (e.g. TMT, E-ELT), 
with ALMA and from the space with JWST will be crucial to make 
an observational breakthrough; one of the main goals of these facil- 



Faint Lyman-Break galaxies 1 5 



ities is to perform spectroscopic surveys of faint sources, reaching 
easily zggo ~ 27. Not only the precise estimation of the redshift 
will be possible, but the validation of the dropout technique at that 
level and information on the total mass (from the stellar velocity 
dispersion), mass distribution and rotation curve, metallicity, dust 
content, galactic wind energetics and stellar population distribu- 
tion will be straightforward to derive from 30-40m spectra of the 
galaxies already discovered in the 8- 10m surveys. The rich phe- 
nomenology revealed by these spectra will give us new insight on 
the formation of galaxies, that will complement the theoretical ef- 
forts made in the next years. 



6 SUMMARY 

We have compared the prediction s of MORGANA+GR ASIL 
( [Monaco et al]|2007l : ISilva et alj[T998l ; iFontanot et ai]|2007h with 
data of Lyman-break galaxies at z ~ 4 — 6. Model predic- 
tions have been produced by simulating deep ACS-HST fields, 
taki ng into account photometric errors as in the GOODS-S sur- 
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vey jCiavalisco et alj[2004l) . and the effect of Lyman-alpha emis- 
sion line. They have been compared to number c ounts of B- V- 
and z- dropouts, scaled to the HUDF sensitivity bv lBouwens et al.l 
( l2007h , colour evolut ion, redshift distribut ions and selection func- 
tion of GOODS-S jVanzella et all llOOgh and GOODS-MUSIC 
jGrazian e t al. 2006^ gal axies, and luminosity functions estimated 
bv lBouwens et al.. (.20071) . 

Model galaxies closely resemble the observed ones in terms of 
sizes, stellar masses, colors, dark matter halo properties and spatial 
distribution but have significantly higher metallicities. However for 
such galaxies, observed at rest-frame UV wavelengths, dust atten- 
uation is as important as uncertain. We find that reasonable tuning 
of GRASIL parameters, namely escape time tcsc of stars from the 
parent molecular cloud and fraction /moi of gas in the molecular 
phase, can result in more than a magnitude difference. The two pa- 
rameters are also degenerate, i.e. decreasing the first has a similar 
effect as increasing the second and vice versa. These parameters are 
also degenerate with model parameters that regulate the amount of 
star formation in high-redshift halos. 

Fixing the parameters it is possible to reproduce the bright 
end of the number counts and luminosity functions, with a mod- 
est overestimate of the number of bright B-dropouts and a possi- 
ble underestimate of bright i-dropouts. We then notice an excess 
of faint (z850 ^ 27) star-forming galaxies, especially when con- 
sidering V-dropouts at 4.5 < z < 5.5. This excess is unlikely 
to be due to sample incomp leteness or to some pe culiar behaviour 
of GRASIL dust attenuation. iFontanot et alj ( l2009h . found a serious 
discrepancy between models of galaxy formation (including MOR- 
GANA) and data for galaxies in the mass range 10^ - lO^"-'^ M0 at 
z = 0. These form too early in the models, are too passive at z < 3 
and host too old stellar populations at low redshift. This clearly im- 
plies an excess of small star-forming galaxies at high redshift that 
forces modelers to suppress star formation at later time, so as to 
recover the correct number of small galaxies at 2 = 0. 

We define a class of galaxies which are over-produced by the 
model. These faint Lyman-break galaxies, forming stars at < 10 
M0 yr~^, have masses that grow in stellar mass from 2 ~ 6, where 
their mass is ~ 3 x 10* — 3 x lO'"* M©, to 2 ~ 4, where their mass 
has grown to ~ 10^ - 10^° Mq. They are hosted in relatively 
massive dark matter halos, with masses in the range 10^^ — lO'^^ 
Mq and circular velocities in excess of 100 km s~^ . These galaxies 
show a stellar mass - gas metallicity relation which is too high 



compared to the (sparse) available evidence. At z = roughly half 
of them are found in small galaxies, the other half having evolved 
into more massive objects. 

We conclude that some feedback mechanism should be at play, 
causing star formation in these galaxies to be suppressed by strong 
gas and metal outflows. These metals would then be important 
polluters of the IGM, and because this pollution takes place be- 
fore the peak of cosmic star formation, constraining this process 
is of great importance also for IGM studies. As originally shown 
bv lDekel & Silkl fl986h . the amount of ejected mass mainly scales 
with the halo circular velocities, and this is true for most semi- 
analytic galaxy formation models. In this way, under very simple 
assumptions, a suppression of these galaxies implies a suppression 
of star formation in those 10^^ Mq halos for which the peak of 
efficiency of star formation is expected at z = 0. Solving this dis- 
crepancy requires either some form of feedback that, surprisingly, 
is most effective at the highest gas densities, or some more exotic 
solution, like higher SN energy in dense environments with low 
metallicity, a top-heavy IMF or a not perfectly cold and collision- 
less dark matter which suppresses such compact halos. Observa- 
tional breakthroughs able to clarify the phisical processes responsi- 
ble for the suppression of these "excess" galaxies may be obtained 
with the future 30-40m optical telescopes, ALMA and JWST; a 
strong theoretical and numerical effort on this topic should be car- 
ried out in the next future to produce predictions that can be tested 
by future observations. 
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